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Abstract 

The Color Glass Condensate picture of the nuclear wave function at 
small- x successfully predicted the suppressed production of high-p^ 
particles at forward rapidities in deuteron-gold collisions at RHIC. This 
triggered more efforts which resulted in theoretical improvements and 
predictions for different observables which will provide further phe- 
nomenological tests. I review recent theoretical developments and dis- 
cuss the resulting predictions. 

1 Saturation and the Color Glass Condensate 

When probing small distances inside a hadron or nucleus with a hard process, one resolves their 
partonic constituents. Increasing the energy of the scattering process at a fixed momentum trans- 
fer allows to probe lower-energy partons, with smaller energy fraction x. As the parton densities 
in the hadronic/nuclear wavef unction grow with decreasing x, they eventually become so large 
that a non-linear (yet weakly-coupled) regime is reached, called saturation, where partons do not 
interact with the probe independently anymore, but rather behave coherently. 

The Color Glass Condensate (CGC) is an effective theory of QCD [1] which aims at de- 
scribing this part of the wave function. Rather than using a standard Fock-state decomposition, 
it is more efficient to describe it with collective degrees of freedom, more adapted to account for 
the collective behavior of the small-re gluons. The CGC approach uses classical color fields: 

\h) = \qqq) + \qqqg) + ■ ■ ■ + \qqqg ■ ■ ■ ggg) + ■ ■ ■ \h) = J DA $ XA [A] \A) . (1) 

The long-lived, large-x partons are represented by a strong color source p^l/gs which is static 
during the lifetime of the short-lived small-x gluons, whose dynamics is described by the color 
field A ~ 1/gs- The arbitrary separation between the field and the source is denoted xa- 

The CGC wavefunction <& XA [A] is the fundamental object of this picture, it is mainly a 
non-perturbative quantity, but the xa evolution can be computed perturbatively. Requiring that 
observables are independent of the choice of xa, a functional renormalization group equation can 
be derived. In the leading-logarithmic approximation which resums powers of as ln(l/xA), the 
JIMWLK equation describes the evolution of |<I> XA [^l]| 2 with xa- The information contained in 
the wavefunction, on gluon number and gluon correlations, can be expressed in terms of n-point 
correlators, probed in scattering processes. These correlators consist of Wilson lines averaged 
with the CGC wavefunction, and resum powers of gsA. 
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Fig. 1: Forward particle production in d+Au collisions at RHIC. The left plot shows the importance of including both 
the large-a; DGLAP evolution of the dilute deuteron and the small-a; CGC evolution of the dense nucleus. The right 
plots shows the excellent description of the spectra shapes, and the K factors needed to obtain the normalization. 



The JIMWLK equation reduces to a hierarchy of equations for the correlators. Most of the 
phenomenology uses a mean-field approximation which significantly simplifies the high-energy 
QCD evolution: it reduces the hierarchy to a single closed non-linear equation for the two-point 
function (S xy ) XA = 1 — (T xy ) XA , the Balitsky-Kovchegov (BK) equation [2,3]- It reads 

d{S xy } x _ a /\ 2 . ._ . ,_ . \ _ (x-y) 



dln(l/x) 



^ J d 2 z M xyz ((S xz ) x (S zy ) x - (S xy ) x ) , M xyz = (x _ ( ^ )2( ^_ y)2 , (2) 



with a = asN c /ir. All the correlators can then be expressed in terms of the solution of this 
equation. Finally, the Fourier transform of the dipole correlator J d 2 xd 2 y e lk '( x-y ) {T xy ) XA / (x— y) 2 
is an (all-twist) unintegrated gluon density. It determines forward particle production [4], while 
more exclusive final states involve more complicated correlators. Solving eq. (|2]) reveals the 
existence of an intrisic momentum scale in the nuclear wave function: the saturation scale Q s {x) 
which characterizes the transition from the dilute regime k > Qs to saturation regime k < Q s . 

One of the most important progress is the recent calculation of the next-to-leading evolu- 
tion equation [5,6]. Concerning how the running coupling should be included, two schemes have 
been proposed by Balitsky (B) and Kovchegov and Weigert (KW). The following substitution 
should be done in formula §2}, with R 2 (x, y, z) given in [6]: 



q s ((x-z) 2 ) 2 a s ((x-z) 2 )a s ((z-y) 2 ) (x-z) • (z-y) a s ((z-y) ; 
(x-z) 2 a s (R 2 (x,y,z)) (x-z) 2 (z-y) 2 (z-y) 2 



a kw N c 
IB 

N c a s ((x-y) 2 ) f 1 / a,((x-z) 2 ) \ 1 / a s ((z-y) 2 ) 

2?r 2 [ xyz+ (x-z) 2 V«s((z-y) 2 ) 7 (z-y) 2 W(x-z) 2 ) 



(3) 
(4) 



At next- to-leading order, there remains a discrepancy between the linear part of the BK equation 
and the BFKL equation. Running coupling corrections to particle production have also been 
investigated [7]. Another important recent theoretical development is the inclusion of Pomeron 
loops in the evolution [8], and the derivation of potential consequences at very high energies 
[9,10]. Concerning phenomenology at present colliders, there was however no significant impact. 




Fig. 2: Angular correlations of Mueller-Navelet jets in the NLL-BFKL framework. Left plots: for standard Tevatron 
(top) and LHC (bottom) kinematics, the renormalization scale uncertainty is indicated (with Q 2 = fcifo). Right plot: 
predictions for CDF (miniplugs) kinematics. The A$ distribution is peaked around A$ = 0, which is indicative of 
jet emissions occuring back-to-back, and the A$ distribution flattens with increasing A 77 or with R deviating from 1. 

2 Forward particle production in pA collisions 

Forward particle production in pA collisions allows to investigate the non linear QCD dynamics 
of high-energy nuclei with a probe well understood in QCD. Indeed, while such processes are 
probing small-momentum partons in the nuclear wavefunction, only high-momentum partons of 
the proton contribute to the scattering (■\fsx p = ke y and ^/sxa = ke~ y ). The dilute hadron con- 
tributes via standard parton distribution functions while the CGC is described by its unintegrated 
gluon distribution. It was not obvious that the CGC picture (Q}, which requires small values of 
xa, would be relevant at present energies. However, it has been the case for many observables 
in the context of HERA [11] and RHIC [12]. One of the most acclaimed successes is the pre- 
diction that the yield of high-p^ particles at forward rapidities in d+Au collisions is suppressed 
compared to A pp collisions, and should decrease when increasing the rapidity. 

In Fig.l the dAu — ► hX pt spectra computed in the CGC approach [13] is compared to 
RHIC data, and the description of the slope is impressive. The need of K factors to describe 
the normalization could be expected since this is a leading-order based calculation. Improving 
the calculation with the next-leading evolution has yet to be done. While the suppression was 
predicted in the CGC approach, other postdictions later offered alternative descriptions. The idea 
is that the value of x probed in the deuteron is so high that large-x effects could be responsible 
for the suppression [14]. This would not happen in pA collisions at the LHC, with a smaller x p . 

While the CGC framework was quite successful in describing single inclusive particle pro- 
duction at forward rapidities, the focus should now shift towards more exclusive observables like 
two-particle production pA—> hxh^X. In particular the correlations in azimuthal angle between 
the produced hadrons should be suppressed compared to pp collisions [15]. By contrast with 
single particle production, in two-particle production the CGC cannot be described only by its 
unintegrated gluon distribution, the so-called -factorization framework is not applicable. This 
means that more tests could be done, probing the CGC structure deeper. The second d+Au run 
at RHIC gives the opportunity to carry out such measurements. 
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Fig. 3: Two-particle production at forward rapidities in pA collisions. The A(j} spectrum is displayed at RHIC (left) 
and LHC (right) energies. When varying pr 2 at fixed j/a, the correlation in azimuthal angle is suppressed as pr 2 gets 
closer to the saturation scale. At the LHC, smaller values of xa are probed, and the azimuthal angle correlation is 
more suppressed as indicated by the vertical axis. The peak is also less pronounced. 



3 Selected predictions for the LHC 

Mueller-Navelet jets [16] in hadron-hadron scattering are two jets produced in each of the for- 
ward directions. In the high-energy regime, in which the jets are separated by a large rapidity 
interval, this process is sensitive to the small-x QCD evolution. An interesting observable is the 
azimuthal decorrelation of the jets as a function of their rapidity separation A77 = 2/1—1/2 and of 
the ratio of their transverse momenta R = £2/^1 [17]. Predictions are shown in Fig. 2, for Teva- 
tron and LHC kinematics, where A<I> = tt — <pi + 4>2 is the relative azimuthal angle between the 
two jets. The curves are obtained in the linear regime, using next-to-leading logarithmic (NLL) 
BFKL evolution. At higher energies, saturation effects will also be relevant [18]. 

Coming back to forward particle production in pA collisions, predictions for the process 
pA — ► hihzX are shown in Fig.3, for RHIC and the LHC [19]. fei, &2 and yi, j/2 are the 
transverse momenta and rapidities of the final state hadrons, and the azimuthal angle spectra are 
displayed. It is obtained that the perturbative back-to-back peak of the azimuthal angle distri- 
bution (which is recovered for very large momenta) is reduced by initial state saturation effects. 
As the momenta decrease closer to the saturation scale (Q s ~ 2 GeV), the angular distribution 
broadens. But at RHIC energies, saturation does not lead to a complete disappearance of the 
back-to-back peak. 

Finally, predictions for the total charged-particle multiplicity in AA collisions at the LHC 
are shown in Fig.4. Two approaches are compared: in the first, kx -factorization is assumed but 
the evolution of the unintegrated gluon densities is accurately obtained from the next-leading 
BK equation [20]; in the second, the x evolution is only parametrized but multiple scatterings 
are correctly taken into account by solving classical Yang-mills equation [21]. While full next- 
leading treatment of both multiple scatterings and small- x evolution is desirable, the numbers 
obtained are similar, which indicates that the uncertainties in both approaches are under control. 
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